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ABSTRACT 

We present the first results of a multi-wavelength survey, incoporating Herschel- 
SPIRE, Spitzer, GALEX and ATCA observations, of a l°xl° field centred on Cen- 
taurus A. As well as detecting the inner lobes of the active galactic nucleus (AGN) 
jet and counterjet, we have found two clouds, bright at sub- mm wavelengths, ~15 
kpc from the centre of Cen A that are co-aligned with the jets. Flux measurements 
at Herschel wavelengths have proved vital in constraining fits to the Spectral Energy 
Distributions (SEDs). The clouds are well fit by a single-temperature, modified black- 
body spectrum (/3 = 2) indicating that we are looking at two cold dust clouds on the 
outskirts of Cen A. The temperature and masses of the clouds are: Tnorth = 12.6l]^'2 
K, Tsouth = 15.ltl-i K; log(Mnorth/M0) = 5.8t^-J, log(Msouth/M0) = 5.6t^-J and 
the gas-dust ratio for both clouds is ~100. The measured values for the northern dust 
cloud are consistent with previous measurements from 750 while the southern cloud 
is a new sub-mm detection. The two dust clouds are located at the termini of the 
partial Hi ring that surrounds Cen A which is also where the gas column density 
peaks. The Herschel survey encompasses the partial H I ring yet we find no evidence 
of dust emission in any other part of the ring. Assuming that the gas-dust ratio is the 
same in the rest of the ring, dust mass upper limits in the H I ring are consistent with 
low column density dust being present but falling below the SPIRE detection limit. 
We have discussed the origin of these clouds and various possible heating mechanisms. 
The observations favour a scenario in which the gas and dust were once part of a late- 
type galaxy which has since merged with Cen A. The dominant heating mechanism 
which adequately explains the observed temperatures in both clouds is heating from 
the evolved stellar population within Cen A. 

1 INTRODUCTION 

1.1 Cen A and its environment 

Centaurus A (NGC 5128) is the dominant, elliptical galaxy 
at the heart of the Centaurus group, 3.6 Mpc away (Chat- 
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topadhyay et aL, 2009). It exhibits a complex variety of gas, 
dust and stellar structures, observable throughout the elec- 
tromagetic spectrum, on scales ranging from sub-pc to kpc 
scales. 

Optical broadband images are dominated by an aging 
stellar population (Soria, 1996), arranged in a system of 
shells (Peng et al. 2002), extending out to ^15 kpc. Bisecting 
the galaxy along its minor axis is a prominent dust lane, first 
observed by John Herschel (1847). The optical shells and the 
presence of a large amount of dust are unusual properties for 
an elliptical galaxy and are indicative of a past merger with 
a late-type galaxy, rich in dust and gas (Baade & Minkowski, 
1954; Malin, 1983). 

Optical spectroscopy of the interstellar medium has re- 
vealed star forming regions and filaments of highly excited 
gas (Blanco et a/., 1975; Peterson et a/., 1975; Morganti et 
al, 1991; Morganti et al, 1992; Rejkuba et al, 2002). The 
high ionisation states have been interpreted as either due to 
a beamed ionising source emanating from the nucleus (Mor- 
ganti et aL, 1991) or shocks being driven into the Interstellar 
Medium (ISM, Bicknell, 1991; Sutherland, Bicknell k Do- 
pita, 1993; Fragile, et al, 2004). 

At radio wavelengths, spectral line observations have 
detected atomic, molecular and ionised gas in the dusty 
disk, while wide- field observations in the near vicinity of 
Cen A have revealed gas in a huge broken ring surround- 
ing the galaxy at a projected distance of ^15 kpc. The ring 
contains substantial amounts of atomic (Shiminovich et aL, 
1994; Oosterloo & Morganti, 2005; Struve, 2010) and molec- 
ular gas (Charmandaris et a/., 2000) and the velocity struc- 
ture varies smoothly over the extent of the ring, consistent 
with rotation about the centre of Cen A. An exception to 
this trend in velocity can be found at the northern terminus, 
where Oosterloo & Morganti (2005) have discovered highly 
turbulent atomic gas. 

Radio continuum images reveal huge lobes extending 
north and south from the galaxy over ^10°. The radio emis- 
sion is dominated by synchotron emission from a jet origi- 
nating in the AGN of Cen A and this is driving plasma at 
relativistic speeds into the surrounding intergalactic medium 
up to distances of 100s kpc. This ongoing interaction results 
in shock heating of the ISM, clearly visible as extended, 
knot-like features in X-ray images (Kraft et a/., 2003, 2007, 
2009; Croston, 2009). Although the northern and southern 
jets exhibit symmetry at large scales, there are crucial differ- 
ences between them. Both have an inner lobe extending to 
^5 kpc and an outer lobe extending to 100s kpc. The differ- 
ence lies in the existence of a northern middle lobe (NML), 
which is connected to and powered by the central engine via 
a large-scale jet (Junkes et a/., 1993; Morganti et a/., 1999; 
Kraft et aL, 2009; Hardcastle et aL, 2009). In contrast, the 
southern jet appears to terminate at the boundary of the 
inner southern lobe where it is driving a strong shock into 
the ISM (Kraft, 2003, 2007; Croston, 2009). 

One aspect of AGN that has come under close scrutiny 
in recent years is their supposed ability to regulate star for- 
mation in the galactic environment by injecting energy into 
the interstellar medium (ISM). This influx of energy pre- 
vents atomic gas from cooling and condensing into molecu- 
lar clouds from which stars would otherwise form (e.g. Silk 
& Rees, 1998). This theoretical mechanism has been invoked 
to try and recover the apparent discrepancy between the ob- 



served population of large, blue galaxies and those predicted 
by cosmological simulations (Benson et aL, 2004; Begelman 
& Ruszkowski, 2005; Binney, 2005). Yet it has also been 
demonstrated that AGN have the ability to create stars via 
the interaction of the jet with gas in the near vicinity of 
the AGN host galaxy (Rees 1989; van Breugel et a/., 2004 
and references therein). In terms of affecting a galaxy's evo- 
lution, these two phenomena seem to be at odds with one 
another and clearly the net effect on the galaxy's star for- 
mation rate (SFR) will depend on the relative efficiencies 
of these two processes. Cen A has been the focus of much 
debate in this respect, with numerous lines of evidence in 
support of jet-induced star formation taking place in Cen 
A (Blanco et a/., 1975; Graham, 1998; Mould et a/., 2000; 
Rejkuba et al, 2002; Morganti & Oosterloo, 2005). 

Observations of Cen A have recently been extended into 
the sub-mm with the detection of a cold dust cloud ~ 15 
kpc north of the nucleus (Stickel et a/., 2004). Using the 
ISO observations at 90, 150 and 200 |J,m Stickel et al. were 
able to model a SED applying a single-temperature modified 
blackbody spectrum (/3 = 2), which is typical of thermally- 
emitting cold dust. The resulting fit implied the dust cloud 
has a very low temperature (13 ^ T ^ 15.5 K) and at the 
distance of Cen A, an inferred dust mass range of 4.5 ^ 
\og{Mdust) ^ 4.9. Unfortunately the peak of the SED for 
very cold dust (< 15 K) occurs longwards of 200 \im. With 
only two detections and an upper limit at 90 ^im, the fit was 
constrained to the rising side of the SED. One of the niches 
that Herschel fills, that we exploit in these observations, 
is that its particular wavelength coverage (70^im- 500^im) 
enables observations of both sides of the peak, thus allowing 
for stricter constraints on the SED fit to the dust thermal 
emission. 

In their discussion of the origin of the cloud, Stickel et 
al. combined their dust mass with the gas mass estimated 
from Schiminovich et al. (1994) and Charmandaris et al. 
(2000) and found that the gas-dust ratio for this cloud is 
approximately 300. This value falls within the range exhib- 
ited by late- type galaxies (Sodroski et a/., 1994; Stickel et 
al., 2000). The ratio of atomic gas to molecular gas is close 
to unity, which is also consistent with values observed in 
spirals (Charmandaris et a/., 2000). This points towards the 
late-type galaxy, with which Cen A is believed to have re- 
cently merged and is responsible for the dust and gas in the 
centre of Cen A, as also being the source of the gas and dust 
observed in the external ring. 

Stickel et al. discuss at length the heating source of 
this cloud. In particular they focus on the ability of stel- 
lar light from Cen A to heat the cloud directly. Temi et al. 
(2003) developed a model for predicting the temperature 
of different types of dust grains at various distances from 
the centre of early- type galaxies. The model accounts not 
only for reprocessed stellar emission, but also for electron- 
grain interactions one would expect in a hot ISM, typical of 
early-type galaxies. One of the main results was that grains 
of different sizes are dominated by different heating mech- 
anisms. Small grains (< 0.1 ^im) are predominantly heated 
by electron- grain collisions in the hot gas, while large grains 
(^1 \im) are heated by the stellar emission of the evolved 
stellar population. 

This model was then used by Stickel et al. to demon- 
strate that the temperature of the dust in the northern cloud 
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is fully consistent with the temperature predicted from the 
Temi et al model, applied to NGC 4636. Since NGC 4636 
is similar in morphology, size and mass to Cen A (Israel, 
1998), they concluded that the FIR emission is mostly due 
to dust heated by the ambient starlight from Cen A without 
the need for heating by other sources such as embedded star 
formation. 

Since their survey field did not encompass the whole 
Hi ring, they were unable to probe the existence of dust 
elsewhere in the gas ring. Also the resolution of ISO was 
not high enough to allow them to comment on the apparent 
alignment between the dust and the AGN jet. The superior 
sensitivity, mapping speed and resolution of Herschel present 
us with the opportunity to map the entire Hi ring in the 
sub-mm to unprecedented levels of sensitivity and detail. 
For more information on Cen A, Israel (1998) provides a 
comprehensive review, while Morganti (2010) complements 
this picture with discoveries that have been made in the 
interim. 

1.2 Very Nearby Galaxy Survey 

The Very Nearby Galaxy Survey (VNGS, P.I.: C. Wilsor£]) 
is a Herschel Guaranteed-Time Key Project focussing on 
twelve galaxies within 25 Mpc and the archetypal starburst 
galaxy, Arp 220. The galaxies show a diverse range of masses 
and properties, from low-mass, late-type galaxies such as 
NGC 2403 to the large elliptical Cen A, with its energetic 
active nucleus. 

The choice of galaxies was also driven by the wealth 
of ancilliary data that are available. The design of the sur- 
vey incoporates both photometry (SPIRE, PACS) and spec- 
trometry (SPIRE-FTS, PACS spectrometer). With such a 
great breadth of multi-wavelength data the survey aims to 
gain a detailed understanding of the processes that regulate 
the ISM, and how these processes vary with the environ- 
ment within different galaxies. The detailed study of these 
resolved galaxies will not only act as a benchmark for stud- 
ies of more distant galaxies, but will also bridge the gap be- 
tween surveys of distant objects and the extensive Galactic 
surveys which have superior physical resolution but by their 
nature are limited to observations of one galaxy. Recent re- 
sults from the VNGS include determining the nature of dust 
heating in M81 (Bendo et al., 2010), attempting to separate 
Galactic cirrus emission from extragalactic dust emission in 
M81 (Davies et a/., 2010) and studying the dust distribu- 
tion and heating mechanisms in and around the starburst 
galaxy, M82 (Roussel, 2010; Panuzzo, 2010). While this pa- 
per concentrates on the external dust around Cen A, a de- 
tailed analysis of the dust within Cen A will be discussed in 
Parkin et al. (2011, in prep). 

In Section [2] we detail the Herschel observations and 
the data reduction techniques that were employed. We also 
give a brief summary of each of the ancilliary datasets that 
were used in the analysis. In Section [3] we compare multi- 
wavelength images of the clouds as seen in SPIRE, Spitzer, 
GALEX and ATCA and show the results of SED fitting to 
the two dust clouds. In Section |4] we discuss the origin of the 
dust clouds and the heating machanisms which are powering 
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the sub-mm emission. Finally in Section[5]we summarise our 
findings. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Sub-mm observations 

The Herschel observations of Cen A consist of a 1° xl° area 
centred on Cen A itself and mapped by both SPIRE (250, 
350 & 500 i^m) and PACS (70 & 160 |^m) using Large Scan- 
map mode. At the time of writing, the PACS data reduction 
pipeline was optimised for the analysis of the central dust 
lane in Cen A, rather than the analysis of the dust clouds 
themselves, hence the regions in the final images surrounding 
the clouds suffered from degraded surface brightness sensi- 
tivity. For these wavelengths MIPS maps (described in Sec- 
tion [212]) proved to have lower surface brightness limits and 
so were used in place of PACS. 

For the SPIRE observations, the telescope was driven 
at a scan rate of 30^^s~^ and the observations repeated five 
times. The SPIRE photometer (Griffin et a/., 2010) data 
were processed up to Level 1 (i.e. calibrated bolometer time- 
lines) with a custom made pipeline adapted from the official 
pipeline 0. This Jython script was run in the Herschel Inter- 
active Processing Environment (Ott, 2010). The purpose of 
the pipeline is to remove all instrumental artefacts such as 
glitches, finite bolometer response time, electronic filtering 
and thermal drift, as well as applying astrometry and flux 
calibration. 

The main difference between our pipeline and the stan- 
dard one is that we did not run the default temperatureDrift- 
Correction and the residual, median baseline subtraction. 
Instead we use a custom method called the BRIght Galaxy 
ADaptive Element (BriGAdE) to remove the temperature 
drift. No further baseline subtraction was necessary to bring 
the bolometer baselines to a common level. 

BriGAdE uses the information from the thermistors in 
each array and directly fits the thermistor data to the en- 
tire bolometer timeline (including data where the spacecraft 
was slewing between scans). If both thermistors exhibited 
instantaneous 'jumps' (an artefact where there is a sudden 
DC offset in the timelines) these are either corrected or the 
comparison switched to the slightly less sensitive Dark Pixels 
of the individual array. This approach is hampered slightly 
by the signal of bright sources in the bolometer timelines. 
To suppress their influence, they are automatically removed 
from the fitting process along with samples affected by other 
artefacts ('jumps' and glitches) in the timelines. When a 
choice of thermistor is available (i.e. for the 250 & 500 \im 
arrays), the one providing the best fit is used to subtract 
a scaled version from the bolometer timelines. This method 
improves the baseline subtraction significantly, especially in 
cases where there are large or rapid temperature variations 
during the observations. 

For the final maps, we used the inbuilt naive mapper of 
HIPE. The final maps have pixel sizes of 6'',8''& 12'' at 250, 
350 & 500 ^Lm respectively and instrument-associated noise 
levels of 5, 5 & 6 mJy/beam. 

2 See Griffin (2008) or Dowell (2010) for a more detailed descrip- 
tion of the ofiical pipeline and a list of the individual modules. 
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2.2 Mid- far-infrared observations 

Cen A was observed with the Spitzer observatory using the 
IRAC and MIPS instruments. We used 24, 70, and 160 ^im 
MIPS raw data from astronomical observational requests 
(AORs) 4940288, 4940544, and 4940800 retrieved from the 
NASA/IPAC Infrared Science Archive to produce the MIPS 
images. The data from AOR 4940288 were taken at the 
medium scan rate and cover regions of ^-^20x50^ that in- 
clude the entire optical disc in each wave band. The data 
from AORs 4940544 and 4940800 were taken at the fast 
scan rate and cover regions of ^6x115^ that only cover the 
jets and the centre of the galaxy. The resulting 160 ^tm maps 
from these fast scan data contains gaps, although the 24 and 
70 \im data completely cover the field of view. 

The data were processed using the MIPS Data Analysis 
Tools (Gordon et aL, 2005) along with additional processing 
steps decribed by Bendo et al (2010) and Bendo et a/. (2011, 
in prep.). The additional processing steps include an addi- 
tional flatfielding step for the 24 \im data; the removal of 
latent images from the 24 |j,m data; background subtraction 
in all three wave bands (including zodiacal light subtrac- 
tion in the 24 \im band and short term drift removal in the 
70 and 160 ^im bands); and a nonlinearity correction ap- 
plied to the 70 ^im data. For the background subtraction to 
be performed correctly, it is important to identify regions 
that are not part of the background; we explicitly identified 
the optical disk and the structures seen in the jets in the 
SPIRE images as regions to be excluded when calculating 
background levels. 

The FWHM of the three MIPS beams are 6 arcsec at 24 
|j,m, 18^^ at 70 |am, and 38^^ at 160^im. The rms background 
levels in the final images are 0.0621, 0.463, and 0.553 MJy 
sr~^ at 24, 70, and 160 ^tm, respectively, and the calibration 
uncertainties are 4% at 24 \im (Engelbracht et a/., 2007), 
10% at 70 M^m (Gordon et a/., 2007), and 12% at 160 ^im 
(Stansberry et al 2007). 

IRAC basic calibrated data frames at 3.6 & 8 ^tm from 
AOR 4939008 were taken straight from the NASA/IPAC 
Infrared Science Archive. The MOPEX software (vl8.7.0) 
was used to interpolate and co-add the individual frames 
into one large mosaic. The 3.6 k 8 [im FWHMs are 1.7" & 
2.0" respectively and the corresponding noise levels in the 
final mosaicked images are 0.025 and 0.06 MJy sr~^. 

As noted in Helou et al. (2004) the 8 \im band contains 
contributions from both stellar emission and PAH emission. 
We adopted their technique to isolate the PAH emission: 
foreground stars were identified by their 3.6-8 ^im colours 
- all unresolved objects with 3.6/8 surface brightness > 4.5 
were masked; the 3.6 & 8 ^im images were then multiplied by 
the aperture correction factors for extended sources (0.944 
and 0.737 at 3.6 & 8 \im respectively); finally the stellar 
contribution was subtracted from the 8 ^im band using the 
formula: 



/,(PAH) = /.(8|am) - 0.232/, (3.6|am) (1) 

The final PAH map was then smoothed to the FWHM of 
the SPIRE 250 ]im beam to enhance the extended emission 
and to allow comparison with the 250 \im map. 



2.3 Radio observations 

Deep 20-cm continuum and 21-cm spectral line observations 
of Cen A taken with the Australian Telescope Compact Ar- 
ray (ATCA) already exist in the literature (20cm contin- 
uum: Morganti et al, 1999; Hi: Oosterloo & Morganti, 2005, 
Struve et a/., 2010). The authors were kind enough to make 
the data available for our use and we only summarise the 
observations here. 

The continuum observations were conducted with the 
ATCA in multiple configurations to maximise the brightness 
temperature sensitivity while retaining high angular resolu- 
tion. The on-source time was 12 hours in each configuration. 
The final size of the cleaned beam in the restored image we 
have used in our analysis is 56^^ x26^^ (PA = 2°) and the rms 
noise is 6.7 mJy/beam. 

The Hi observations were also conducted using the 
ATCA using a combination of two configurations. Due to 
the strength of the continuum source in Cen A, they were 
restricted in the use of short baselines, which imposed a 
high-pass filter on the final reduced images. This resulted in 
all emission on scales larger than 7' being filtered out. This 
however does not affect our analysis. The size and shape of 
the restored beam in the final datacube was 29^^x19.5^^ (PA 
= -4.7°). A moment map was produced from the datacube 
with a 3-cr detection limit of l.OxlO^^cm"^ over a velocity 
width of 13.2km s~^. For further details the reader is re- 
ferred to Morganti et al (1999) and Oosterloo & Morganti 
(2005). 

2.4 Ultraviolet observations 

Cen A was observed by GALEX (Martin et a/., 2005) as part 
of the Nearby Galaxy Survey (Gil de Paz, 2004). The ob- 
servations were carried out simultaneously in far- ultraviolet 
{FUV: Xeff = 1539 A, AA = 442 A) and near- ultraviolet 
{NUV: Xeff = 2316 A, AA = 1060 A). The images used 
here were taken from the GALEX GR6 data release and 
have total exposure times of 20101 and 30460 s in FUV and 
NUV. 

The final data products consist of circular images ^-^1.2° 
in diameter with a resolution of 4.5'' (FUV) and 6'' (NUV). 
The surface brightness limits for the FUV and NUV bands 
are, 30.6 and 29.5 mag arcsec"^ respectively. The reader is 
referred to Morrissey et al. (2007) for details of the instru- 
ment, the obervations and data reduction. 



3 RESULTS 

3.1 SPIRE images 

Fig. [T] shows the entire observed field at 250, 350 & 500 
\im. The prominent dust lane that seems to cut the galaxy 
in two at optical wavelengths is clearly seen in emission at 
sub-mm wavelengths in the centre of the field. The inner 
radio lobes are clearly visible as extensions perpendicular to 
the dust disk, particularly in the 350 & 500 \im bands. As 
stated in the introduction, cold dust peaks longward of 200 
\im. This means that the 350 & 500 \im bands are situated 
well on the Ray leigh- Jeans side of the peak of the black- 
body curve. Regions of thermal emission are characterised 
by decreasing flux levels at increasing wavelengths. From 
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Figure 1. Herschel-SPIKE maps of Centaurus A with insets showing the two dust clouds, top: 500 |J.m, bottom: 350 |J.m, overpage: 250 
|xm. The dusty disk of Cen A dominates the emission at the centre of the image. Synchrotron emission from the inner lobes of the AGN 
jet are visible as extensions northeast and southwest from the centre of the galaxy. All colour bars are in units of Jy/beam. 
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Figure 1 - continued Herschel-SPIRE maps of Centaurus A: 250 |xm. 



comparing the lobes at the different SPIRE wavelengths, it 
is clear that the they exhibit the opposite emission trend, 
suggesting that another emission process must be contribut- 
ing to these longer wavelengths. Synchotron emission from 
relativistic electrons spiralling around magnetic field lines 
in the jet is the most obvious candidate. Synchrotron has a 
well characterised spectral slope, a, such that S oc and 
thus the longer wavelength Herschel maps can have signifi- 
cant contributions from synchrotron emission. Synchrotron 
will be discussed in the section [331 since its presence will ad- 
versely affect the SED fitting routine, which only constrains 
a thermal spectrum. 

Following the line traced by the jets, we see the two 
clouds: -15' NE from the centre of Cen A (13^26^ 18.876^ 
-42°49'32'') there is a knot of emission about 4' across 
and -12' SW from the centre of Cen A, (13^24^34.976", 
-43°09'02'') there is another cloud -2' across. Unlike the 
emission from the inner lobes, the clouds are brighter at 
shorter wavelengths, suggesting that the sources have a sub- 
stantial thermal component. The northern cloud is undoubt- 
edly the dust cloud that was detected by Stickel et al (2004), 
but the southern cloud, to our knowledge, is a new detection. 

The morphology of the clouds is very different; the 
southern cloud is fairly smooth with an elliptical shape, 
whereas the northern cloud is more flocculent and appears 
to be two smaller clouds separated by a dust-free lane that 
runs almost perfectly N-S. 

The Cen A survey field is large enough that the zero 
level can be affected by the presence of Galactic cirrus, which 
introduces gradients across the image. Naturally this com- 
plicates our estimates of the background in each region. Un- 



der the assumption that the emission is optically thin, the 
southern cloud could be contaminated with more substan- 
tial amounts of Galactic cirrus emission than the northern 
cloud. This explains why the southern source itself appears 
brighter, but as we shall see in Section 13.61 the integrated 
flux of both clouds is actually fairly similar in each band. 



3.2 Spitzer images 

At other wavelengths the clouds show some interesting fea- 
tures. Figs. [2H5] illustrate how the mid-infrared (MIR) and 
far-infrared (FIR) structures compare with those detected in 
the sub-mm. Each image is a 5'-10' region immediately sur- 
rounding each cloud and shows the 250 ^im contours overliad 
on the images at each individual wavelength. Since 160 ^tm is 
close to the peak of the SED for cold dust, one would expect 
that the morphology exhibited at 250 ^im should be similar 
to that at 160 ^im and this is what we find. In Fig. [2] the 
clouds closely trace the 250 ^im emission, although the loca- 
tion of the northern cloud places it close to the edge of the 
area covered by the MIPS medium scan rate image, hence 
the area around the northern cloud contains some artefacts 
related to edge effects at the end of the medium scan rate 
maps and blank pixels from regions not covered in the fast 
scan maps. The southern cloud is a strong detection and the 
extended emission seen at 250 ^im is also seen at 160 ^im. 

If the dust is cold, the spectrum drops off rapidly at 
wavelengths shorter than —160 \im. The expected 70 ^im 
emission is a few orders of magnitude lower than that at 
160 ^im so a lack of any detection at 70 ^im would indicate 
a cold object. While it seems there is a peak in Fig. [3] at 
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1 3^25^00" 24™50" 24"M0" 24"'30" 24™20" 



Eight Ascension (JSOOO) Ri^M Ascension aSOOO) 

Figure 2. Clockwise from left: Spitzer 160 |J.m map of Cen A region. The colourmap has been stretched to enhance diffuse, extended 
emission. There are a number of visible artefacts, but both dust clouds are visible at the centre of the highlighted regions; the region 
around the northern cloud with 250 |xm contours overlaid at 5, 15, 30 & 50 mJy/beam; the region around the southern cloud with 250 
|Ltm contours overlaid at 15, 30 &; 50 mJy/beam. The 160 |Lim emission clearly traces the 250 jiim emission. 




1 3^25^00" 24'"50" 24'"40" 24'^30" 24=^20" 



Right Ascension (JSOOO) Ri^M Ascension (JSOOO) 

Figure 3. Clockwise from left: Spitzer 70 |Lim map of the Cen A region. The colourmap has been stretched to enhance diffuse, extended 
emission. There are significant stripes present which are artefacts from the observing system that could not be removed during the data 
reduction process. One of the stripes extends to the position of the southern dust cloud, making flux measurements highly uncertain; 
the region around the northern cloud with 250 jam contours overlaid at 5, 15, 30 & 50 mJy/beam; the region around the southern cloud 
with 250 |J.m contours overlaid at 15, 30 &; 50 mJy/beam. No 70 |J.m emission can be unambiguously associated with either cloud. 



© 2007 RAS, MNRAS 000, [T]-?? 



8 Auld et al 




13^25^00" 24'"50" 24™40" 24™30" 24^20" 



Eight Ascension (JSOOO) Ri^M Ascension aSOOO) 

Figure 4. Clockwise from leftSpitzer 24 |j.m map of the Cen A region. The colourmap has been stretched to highhght diffuse extended 
emission. The northern inner lobe of the jet is just visible as are the dust clouds; sub-image of the northern dust cloud with SPIRE 250 
|xm contours overlaid at 5, 15, 30 & 50 mJy/beam; sub-image of the southern dust cloud with 250 |xm contours at 15, 30 &; 50 mJy/beam. 
Both the dust clouds exhibit associated peaks at 24 |Ltm. 




Right Ascension (J2000) 



Figure 5. left: Spitzer 8 |J.m map of the region south-west of Cen A. The colourmap has been stretched to enhance diffuse, extended 
emission. The stellar emission from Cen A dominates the north-eastern corner. The southern dust cloud is visible near the centre of the 
image; stellar-subtracted 8 fxm map of the southern dust cloud with SPIRE 250 |J.m contours overlaid. The map has been smoothed to 
the 250 |xm beam and the 250 |xm contours are: 15, 30 & 50 mJy/beam. Blank regions indicate foreground stars which have been masked. 
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Figure 6. Multi wavelength observations of the immediate surroundings of Cen A. SPIRE 250 |J.m false colour image with 20cm continuum 
emisssion (white contours) and Hi column density (black contours). The image colour scale has been stretched to enhance low surface 
brightness features, contour levels are as follows: 20cm - 12.5 x 1, 2, 3, 4, 5, 6, 7, 8 &; 9 mJy/beam; Hi- 1, 5, 10, 15 & 20 xlO^o cm~2. 
The image is centred on Cen A and highlights the regions of the H I ring. 
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Figure 7. Clockwise from left: GALEX NUV map of the Cen A region. The map has been smoothed to the SPIRE 250 |xm beam 
size and the colour scale has been stretched to enhance low-surface brightness, extended features. There is a lot of extended emission 
corresponding to the stellar population of Cen A. A jet-like feature can be seen extending far beyond the stellar radius to the north of 
Cen A; sub-image of the northern dust cloud with SPIRE 250 |J.m contours overlaid at 5, 15, 30 Sz 50 mJy/beam. The region exhibits 
a large amount of UV coincident with the dust emission; sub-image of the southern dust cloud with 250 |J.m contours at 15, 30 &; 50 
mJy/beam. This dust cloud has no detectable UV emission associated with it. 
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Figure 8. Clockwise from left GALEX FUV map of the Cen A region. The map has been smoothed to the SPIRE 250 |J.m beam size 
and the colour scale has been stretched to enhance low-surface brightness, extended features. As well as the extended emission from Cen 
A, Galactic cirrus can be seen in the eastern half of the map. The jet-like feature noticed in the NUV is also visible here; sub-image of 
the northern dust cloud with SPIRE 250 |J.m contours overlaid at 5, 15, 30 &; 50 mJy/beam. The region exhibits a large amount of UV 
coincident with the dust emission; sub-image of the southern dust cloud with 250 |xm contours at 15, 30 & 50 mJy/beam. This dust 
cloud has no detectable UV emission associated with it. 



the location of the southern cloud, it happens to be on top 
of a streak in the image which is a data reduction artefact. 
Because of this, we consider the southern detection spurious 
and have only used data from the 70 maps to assess 
upper limits. Using the rms in the vicinity of the clouds 
leads to an upper limit of ^200 mJy. 

One of the most interesting results to come out of the 
Spitzer Nearby Galaxy Survey was the tight correlation be- 
tween the 24 \im emission and the atomic hydrogen spec- 
tral line emission from Paschen-a, Pa (Calzetti et al, 2005, 
2007; Prescott et a/., 2007; Zhu et a/., 2008, Kennicutt et a/., 
2009). This suggests that the 24 \im band is a particularly 
good tracer of star forming regions (although see Bendo et 
aL, 2010 for alternative sources of 24 \im emission). In Fig. 
H]we see that both of the clouds show 24 \im emission coin- 
cident with the peaks in the 250 \im emission, suggesting a 
common source for both bands. 

Emission from galaxies in the 8 ^im band of Spitzer is 
dominated by broad line emission from polycyclic aromatic 
hydrocarbons (PAHs). There is also a correlation between 
the 8 \im PAH emission (corrected for stellar emission) and 
Pa, however the correlation is not as tight as for 24 \im. The 
relationhip is non-linear and more dependent on metallicity 
than the 24 \im band, so it is not as reliable to associate 8 
\xm emission with star formation (e.g. Calzetti et al, 2005; 
Bendo et aL, 2008). Unfortunately the northern cloud falls 
outside the survey field in the 8 ^im band, but the southern 
field shows some point-like sources embedded in a diffuse, 
extended cloud. When smoothed to the 250 \im beam, the 
point sources blend together so that the emission appears 
as a cloud which closely traces the 250 |J,m contours, again 
suggesting a common origin. 



3.3 Radio images 

Fig.[6]shows the location of the northern jet (white contours) 
and the H i ring (black contours) overlaid on the 250 \im false 
colour image. The systemic velocity of the gas is beyond the 
extent of Galactic H i, which indicates that the gas is part of 
the Cen A system. The dust emission correlates well with the 
H I contours, suggesting that the dust is associated with the 
Hi in Cen A as opposed to being a coincidental projection 
of Galactic cirrus (c.f. M81: Davies et aL, 2010). 

Previous comparisons of the Hi velocity dispersion of 
the southern tip of the northern cloud match the velocity 
range of ionised gas just to the east of the cloud (Ooster- 
loo & Morganti, 2005). This kinematic agreement combined 
with the spatial coincidence is a good indicator that an in- 
teraction is taking place between the jet and the Hi cloud. 
The dust is also coincident with the jet, at least in projec- 
tion. If the dust is well mixed with the gas, as we'd expect, 
it would imply that the dust is also involved in the jet-cloud 
interaction. It is also interesting to note that the western 
part of the H l ring shows no corresponding enhancement in 
the sub-mm. Possible explanations for this are discussed in 
Section HXn 



3.4 Ultraviolet images 

Both of the GALEX images (Figs. El & [8|) have been 
smoothed to the resolution of the SPIRE 250 [im beam to 
help bring out some of the more diffuse UV structure. UV 
emission from the stellar population can be clearly seen in 
both images. A large amount of emission is present in the 
eastern half of Fig. [8] which is mostly likely scattered light 
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from Galactic cirrus. Also striking is the jet-like feature ex- 
tending from Cen A to the north, far beyond the optical 
extent of the galaxy. Focusing on the clouds themselves re- 
veals a stark difference between them. The northern cloud 
clearly has UV emission coincident with it, but the emis- 
sion is mostly restricted to the edges of the cloud. It may 
be a chance alignment, but is also what one would expect 
if the UV emission is embedded within, and obscured by, 
the dust cloud. The southern cloud, however, shows no sign 
of UV emission in either of the GALEX bands. It is possi- 
ble that there is no UV originating from the southern dust 
cloud which would indicate the lack of any star formation in 
the vicinity. Since UV emission is highly attenuated by dust 
grains, this could also indicate that UV is present but has 
been obscured by high column density gas and dust. In the 
next section, we combine the SPIRE, Spitzer and radio data 
to assess whether the dust cloud emission is purely thermal 
or contaminated by synchrotron. 

3.5 Synchrotron Contamination 

Since the inner lobe is particularly bright in the 500 and 
350 \im image, we must consider whether or not the fluxes 
are contaminated by synchrotron radiation from the jet. As 
mentioned in Section [3T] synchrotron radiation has a power- 
law emission spectrum. It is therefore possible to estimate 
the amount of contamination at the SPIRE wavelengths by 
measuring the clouds' fluxes from lower frequencies (e.g. 
from radio maps) and extrapolating to the SPIRE wave- 
lengths. 

Ideally the spectral index would be measured from a 
combination of radio maps, but the map of Morganti et al. 
(1999) is the only radio map with sufficient spatial resolution 
that extends to the distance of the clouds. In this case we 
are reduced to assuming a value for the spectral slope. The 
most extreme case would be to assume that the spectral 
index is as flat as possible, i.e. a — 0.5 (Begelman, Blanford 
& Rees, 1984). We measure a radio flux inside the aperture 
used for the SED fitting of 270±60 mJy at 1.4 GHz. Even 
with the flat spectral index this only produces a flux of a 
few mJy at the 500 |J,m band. Hence, there is no significant 
contamination by the jet in the SED of the northern cloud. 

To our knowledge no radio maps of sufficient resolution 
and sensitivity exist for the region covering the southern 
cloud. Since there is no equivalent to the Northern Middle 
lobe connecting the southern inner lobe to the giant outer 
lobe, the amount of synchrotron must be even lower than in 
the northern cloud. Having rejected synchrotron as a possi- 
ble contaminant we model the fluxes under the assumption 
that the emission is thermal in origin. 

3.6 SED fitting 

Stickel et al (2004) fit the FIR SED for the northern cloud 
using a modified blackbody spectrum, but their analysis was 
limited to two points on the Wein side of the spectrum (150 
& 200 ^im) plus another upper limit (90 \im). With so few 
data points the peak was poorly constrained. With the com- 
bination of MIPS and SPIRE data, we now have 5 points 
straddling the peak, which should enable us to fit a modified 
blackbody in the same manner, but with a higher degree of 
assurance. 



Table 1. FIR-sub-mm flux measurements for the two dust clouds. 
Uncertainties are shown in brackets and are 20% for SPIRE fluxes 
and 30% for Spitzer measurements. 



Cloud S^QQ S'350 *S'250 SiQQ S70 

Jy Jy Jy Jy Jy 



North 0.33 (0.07) 0.8 (0.2) 0.8 (0.2) 0.7 (0.1) < 0.2 
South 0.37 (0.07) 0.9 (0.2) 1.0 (0.2) 1.5 (0.5) < 0.2 



We adopted the SED fitting routine as used in Smith 
et al. (2010), in which the uncertainties in the fitted dust 
mass and temperature are derived from several thousand 
Monte Carlo simulations. The benefit of the simulations is 
the determination of the error distributions for the mass 
and temperature which are not necessarily Gaussian, thus 
lending a higher degree of robustness to the quoted uncer- 
tainties. The fitting routine allows the user to specify the 
emissivity index, /3 as a free parameter, but we found that 
the simulations did not converge, so /S was fixed at 2. The 
dust mass-opacity coefficient used was 0.192 m^ kg-^ at 350 
micron (Draine & Li, 2001). 

Before fitting the SED, the images were smoothed to the 
resolution of the SPIRE 500 \im beam. An aperture was de- 
fined for each cloud based on the highest signal-noise map 
(250 ^Lm) and these apertures were used across all bands. 
The lack of coverage in the 160 ^tm image in the north lim- 
ited our ability to measure flux over the whole cloud, but 
only a small section of the eastern-most part of the cloud 
was lost and it is unlikely to make a significant difference to 
the final result. 

The measured fluxes are shown in Tabled Quoted un- 
certainties are larger than the nominal calibration uncertain- 
ties for each instrument, but given the amount of foreground 
contamination against which the fluxes were measured, we 
feel these larger error bars are justified. It was apparent that 
even after smoothing, the quality of the 70 \im images was 
only good enough to provide upper limits, but this still gave 
stricter constraints than omitting the 70 \im data altogether. 

The SEDs for the dust clouds are presented in Fig. [9l 
The SED for the northern cloud compares well with the SED 
of Stickel et al (2004), and we derive a similar dust temper- 
ature: T = 12.6tl:2 K. We recover more mass than Stickel 
et a/., log(Mdust/M0) = 5.81q;2, most likely because of the 
increased wavelength coverage provided by Herschel. For the 
southern cloud we measure a dust mass and temperature of 
log(Mdust/M0) = 5.6l°;^, T = 15. ll};^ K. Combining the 
dust masses with the gas masses measured by Charmandaris 
et al. (2000), we derive gas-dust ratios of 121 & 110 which 
are typical of late- type, star- forming galaxies. Having es- 
tablished a thermal origin for the sub-mm emission in these 
clouds, we now discuss their origin and heating mechanisms. 
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Figure 9. Modified blackbody, single-temperature fits to the 
FIR-sub-mm data assuming a /3 = 2 emissivity. Data points from 
this analysis are shown as red dots and the data points of Stickel 
et al. are shown as blue squares. Arrows indicate upper limits. 
Top: northern cloud. Bottom: southern cloud. 



4 DISCUSSION 

4.1 Origin of the clouds 

^.1.1 Galaxy merger 

A comparison of the distribution of the atomic gas and the 
dust (Fig. [G]) reveals that the only detected dust emission, 
that can be unambiguously associated with the H i ring, is at 
the position of the two clouds. The lack of sub-mm emission 
in the rest of the gas ring has three possible explanations: 
either there is no dust in other parts of the ring, the col- 
umn density of dust falls below the SPIRE detection limit, 
or the dust is present but too cold to emit at the SPIRE 
wavelengths. 

If the galaxy that merged with Cen A was a late- type, 
it would be expected that the dust is well mixed with the 
gas in the Hi ring. It is unlikely therefore, that there is no 
dust anywhere else in the orbitting gas. This then leaves the 
question of whether the dust is too diffuse or too cold. 

Assuming that the dust in the HI cloud to the west is at 
the same temperature as the two detected clouds, the SPIRE 
detection limit can be used to estimate the limiting column 
density of dust within it. Apertures were used to define the 
local rms in the vicinity of the cloud and this was adopted 
as the noise level. The source region was defined from the 
Hi map of Struve (2010), using an Hi column density of 
10^^ cm~^. The 500 |J,m map was the most sensitive in the 
region of the western cloud. Despite having a slightly larger 
instrumental noise value, this wavelength is less sensitive 
to Galactic cirrus which dominates the background in this 
region. The rms pixel value in the vicinity of the Hi cloud 
was 10.1 mJy/beam. This gave a 3-cr flux upper limit of 



66.6 mJy over the extent of the cloud. For a 14K cloud 
with /3 = 2, this resulted in an upper limit for the dust 
mass of Mdust < 9.4 x lO^M©, or a corresponding gas-dust 
ratio of Mgas/ Mdust > 63. This is not a particularly strong 
lower limit for the gas-dust ratio, and suggests that the dust 
is simply too diffuse to be detected by SPIRE. Without a 
strong upper limit on the gas-dust ratio there is little use in 
exploring the impact of a significant amount of colder dust 
in the western Hi cloud. 

These findings are consistent with having dust well 
mixed throughout the gas ring and the gas-dust ratio is 
typical of late-type galaxies as is the ratio of the atomic 
to molecular gas. This supports the scenario that the dust 
and gas originated in the galaxy that merged with Cen A to 
form the dusty disk at the centre of Cen A. 

The northern cloud is coincident with the high- velocity 
H I gas which has been linked to the jet-induced star form- 
ing region nearby. Opportunities to explore jet-induced star 
formation in the kind of detail that Cen A offers are rare, so 
it is interesting to explore the possibility that this local star 
formation is responsible for the existence of the northern 
dust cloud. 



4' 1-2 Jet-induced Star Formation 

If the dust originated from the stars in the vicinity of the 
clouds, we can use a simple argument to determine whether 
or not the SFRs are sufficiently high to have formed the mass 
of dust we observe. Stardust is formed in the cool, stellar 
winds of low- intermediate mass stars (LIMS) on timescales 
of 0.5 — 1 Gyr and in supernovae (SNe) on timescales of 
tens of Myr. In the Milky Way, the major dust source is 
presumed to be the former, with LIMS injecting dust at a 
rate of -2 xlO^^Moyr"^ (e.g. Whittet 2003). Recent FIR 
observations of Galactic and extragalactic SN remnants with 
Spitzer and Herschel (e.g. Rho et al. 2008; Otsuka et al. 2010; 
Barlow et al. 2010) suggest SNe contribute a dust mass-loss 
rate of < 1 x lO^^Moyr"^ (corresponding to < 0.1 M© 
of dust formed per SNe). However, as SN shocks are also 
responsible for destroying dust in the ISM via sputtering, it 
is not yet clear how much of the dust observed in remnants 
survive their journey into the ISM. 

We can estimate the stellar dust injection rate needed 
to produce the observed dust mass in the Cen A clouds by 
comparing the rate at which dust grains are destroyed in the 
ISM. The timescale for dust destruction can be estimated 
using (McKee 1989): 



^des — 



Ma 



RsNMd.s OmSFR 



(2) 



where Md is the dust mass, Rsn is the SN rate, Md,s is 
the mass of dust shocked and destroyed by one SN blast 
wave (^3 M© for shocks with speeds v > lOOkms"^ - 
Dwek & Cherchneff 2011). Using Eq. [H dust is destroyed 
on timescales of Tdes ^400 Myr in the Milky Way (Jones et 
al. 1994; Tielens 1998). The estimate requires knowledge of 
the SFR, and a number of prescriptions exist to estimate the 
SFR for the northern cloud based on its luminosity at 24 \im, 
NUV and FUV. For completeness the SFRs for the south- 
ern cloud are also calculated, but the reader is reminded 
that there is no evidence of a connection between the south- 
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ern jet and this dust cloud so in this case the assumed star 
formation is not jet-induced. 

Using the prescription in Relaho et al. (2007) as pre- 
sented by Calzetti et al. (2010) the northern cloud exhibits 
5*24 = 9 mJy, which corresponds to a SFR24 = 1.4 x 
lO~^M0yr~^. For the southern cloud we measure S24 = 7 
mJy. This corresponds to a SFR24 = 1.2 x lO"^M0yr~\ 

Iglesias-Paramo et al (2006) provide recipes for cal- 
culating SFRs in the UV bands of GALEX. The UV lu- 
minosities must first be corrected for Galactic extinction. 
Galactic extinction magnitudes, E{B — V), were provided 
by the GALEX, MAST website using objects identified close 
to (< 30^^) the cloud positions within the Gen A field. For 
the northern cloud: E{B — V) — 0.102 and for the southern 
E(B — V) = 0.138. The conversion from Ab to Auv was cal- 
culated using the method described in Wyder et al. (2007). 
The attenuation values are: Auv — 0.840 for the northern 
cloud and Auv = 1.055 for the southern cloud. The mea- 
sured NUV and FUV magnitudes were then corrected by 
these values. 

For the southern cloud, the 5-cr limits based on the er- 
ror in the background estimates, corrected for Galactic ex- 
tinction, were used to provide upper limits. The extinction- 
corrected luminosities were then used to calculate the SFRs. 
For the northern cloud we derive a SFRfc/v = 3.8 x 
lO"^M0yr-^ and SFRivc/v = 2.4 x lO-^Moyr"^ for the 
southern cloud SFRfc/v < 2x lO^^Moyr"^ and SFRatc/v < 
5 X lO~^M0yr~^. All the luminosity measurements and as- 
sociated SFRs are shown in Table [2] For reference, the SFR 
derived from Mould et al. (2000) is also shown. 

The SFR values for the northern cloud as judged from 
four independent tracers differ by factors of a few. The 
SFRs broadly agree, given the different caveats and scat- 
ter that are associated with each SFR estimator. The low 
SFR for the northern cloud implies dust in this environ- 
ment could survive for 2 Gyr. The dust injection rate from 
stars for this cloud would therefore need to be Md/rdes ^ 
10~^ M0 yr~^ to produce the total dust mass. From the low 
SFR lO~^M0yr~^), we would expect a dust injection 
rate from LIMS and SNe (scaled from the Milky Way pa- 
rameters and SFR) of 10~^ M0 yr~^; this is an order of mag- 
nitude lower than required. At this rate, it would take more 
than 30 Gyrs to build up the 4 x 10^ M0 of dust observed 
by Herschel with a stellar mass- loss source of dust. 

The dust clouds detected here cannot be explained by 
replenishment of the dust via stellar mass loss; the destruc- 
tion timescales would need to be more than an order of mag- 
nitude longer than currently estimated and the induced star 
formation would need to be replenishing dust for an unfeasi- 
bly long time. This reinforces the notion that the dust most 
probably originated in the galaxy that merged with Gen A 
and was re-distributed along with the gas as the galaxy was 
torn apart. 

Simulations of merger events have demonstrated that 
the gas and dust are preferentially stripped, so the lack of 
a stellar enhancement in the H i ring is not unusal. It does, 
however, raise the question of what possible heating source 
could be responsible for heating the clouds to produce their 
sub-mm emission. 
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Table 2. IR and UV star formation rates for the dust clouds. 



Cloud North South 



S24 (Jy) 


0.009 


0.007 


L24 (xl039erg s-i) 


1.69 


1.36 


SFR24(M0yr-i) 


0.00143 


0.0012 


"^NUV 


15.3 


> 20.96 


m.FUV 


16.1 


> 19.06 


E{B - V) 


0.102 


0.128 


l^NUV (xl040ergs-i) 


5.92 


< 0.06 


Lfuv (xlO^Oergs-i) 


4.21 


< 0.15 


SFRatc/v (M0yr-i) 


0.0024 


< 5 x 10- 


SFRfuv (M0yr-i) 


0.0038 


< 2 x 10- 


SFRmoo;! (M0yr-i) 


0.0012 





^ Mould et a/. (2000) based on a the maximum age of detected 
OB associations and total number of detected OB associations 
assuming a Salpeter mass function 



4.2 Dust heating mechanisms 

The dust clouds reside in a complex environment. They ex- 
ist on the very outskirts of the stellar extent of Cen A and 
heating by the stellar population will have some effect. The 
region around Cen A is also rich in X-ray emission and the 
northern cloud appears to lie in the projected path of the 
AGN jet and also close to the beam of ionising radiation 
discovered by Morganti et al. (1991). We now consider each 
of these heating mechanisms and discuss their relative con- 
tribution. 



4-2.1 Heating by aging stars 

The model of Temi et al. (2003) predicts temperatures of 
dust grains of different sizes as a function of radius for early 
type galaxies such as Cen A. The model environment con- 
sisted of X-rays from a hot (> 10^ K) ISM and starlight 
from the evolved stellar population. From Fig. 2 of Temi et 
al. (2003) the expected dust temperature at the distance of 
the dust clouds is 14-16 K for grains of size ^ l\im, which 
is entirely consistent with our measurements and also the 
previous measurement of Stickel et al. (2004). For smaller 
grains, the temperature profiles are too high to account for 
the low temperatures we see, so the clouds must be domi- 
nated by large grains if the Temi model is representative of 
the grain population. The temperature profiles from grains 
of this size are dominated by heating from the aging stars, 
with electron-heating becoming more important for smaller 
grains. From this we infer that the clouds are most likely 
heated by the old stellar population. 

It is curious that in such a complicated environment, a 
single heating mechanism can account for the grain heating. 
We now examine the alternative possible mechanisms and 
attempt to explain why they are not credible. 
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4-2.2 Beamed radiation 

lonisation of the material close to the dust clouds in Cen A 
by a beamed nuclear source was first proposed by Morganti 
et al. (1991, 1992). Through extensive modelling incorpo- 
rating multi-cloud photoionisation, they were able to reject 
other scenarios such as heating from a diffuse X-ray halo or 
cosmic ray heating in favour of the beamed, ionising pho- 
tons, emanating from a source embedded deep within the 
obscured nucleus of Cen A. We have used their estimate of 
ionising photon flux to calculate the temperature of a dust 
grain at the distance of the clouds. 

In the absence of attenuation, the total energy incident 
on a grain is simply given by: 

Qheat — QhVHa^grain (3) 

where q is the ionising photon flux (photons s~^ sr""*^), huHa 
is taken as an estimate of the minimum energy of the ionising 
photons and ^grain is the solid angle of a single grain (^3.3 
xlO~^^ sr) at a distance of the dust clouds (~15 kpc). 

This incident power is then absorbed completely by the 
grain, and radiated away: 

Qcool = ^7Vax<^SBTd{Qabs) (4) 

where ax is the cross-sectional area of a grain (a spher- 
ical grain diameter of 1 ^tm is assumed) and Qats is the ab- 
sorption co-efficient and we use the approximation Qabs ~ 
1.35 X IQ-^T^ax following Temi et al. (2003) and Draine & 
Lee (1984). 

By conservation of energy, Eqns. dS]) & dH) must balance. 
So if one knows it is possible to calculate T. Morganti et 
al. (2001) calculate a value of ^ = 10^^ photons s~^ sr~^ to 
account for the Ha luminosity in the filaments. With this 
incident photon flux, it is possible to raise the grains to 
67 K. Since the dust temperature is such a weak function of 
the incident power, one would only require ^0.01% of power 
from the beamed flux in order to raise the dust temperature 
to ^14 K. While this is plausible, the beamed flux would 
have to have a sufficient opening angle to illuminate both 
dust clouds entirely. The northern cloud is roughly 5 kpc 
across and at the distance of the cloud (15 kpc) this would 
require that the photon beam have an opening angle of at 
least ^20°. This is uncomfortably large for a beamed source 
and we note that this minimum estimate does not engulf the 
optical filaments to the east. Also there is no known evidence 
of optical filaments to the south of Cen A, suggesting that 
the ionising beam can't be interacting with the southern 
cloud at all. 

4.2.3 X-ray heating 

X-ray observations of Cen A have revealed that the galaxy 
has an X-ray jet and exists inside a diffuse X-ray halo (Kraft 
et al, 2009; Turner et al, 1997; Feigelson et a/., 1981). The 
model of Temi et al. (2003) takes into account the effects of a 
diffuse X-ray medium in which the dust grains are situated, 
and has been found to impact only small (< 0.1 ^tm) grains. 
It is possible that the presence of discrete X-ray sources 
might have a more significant effect. 

In a recent study, Kraft et al. (2009) reported the dis- 
covery of several bright X-ray knots on the eastern edge of 



the northern middle lobe, which places them to the north 
and east of the northern dust cloud. The nearest X-ray knot 
to the northern cloud (N4 in the notation of Kraft et al.) 
has a luminosity of only ^ 5 x lO^L© and therfore cannot 
be a potential heating source. 

A sufficiently bright low mass X-ray binary (LMSB) 
would also have the potential to heat the cloud. From Figs. 2 
& 3 in Kraft et al. there are two compact X-ray sources close 
to the location of the dust clouds. Using the XMM-Newton 
archive the two sources were located in the original data and 
found to have luminosities of Lx (north) = 2.9 x 10^^ erg s~^ 
and Lx (south) = 3.6 x 10^^ erg s~^. The luminosities are 
orders of magnitude below that required to power the dust 
clouds, so we conclude that X-ray emission either by diffuse 
hot gas or compact sources is not strong enough to account 
for the observed luminosity in the sub-mm. 

4.2.4 Jet mechanical heating 

The final situation to consider is a direct interaction between 
the radio jet and the dust cloud. In their comparison of the 
jet power to the X-ray power of the bright knots, Kraft et 
a/. (2009) estimate that the jet itself has a mechanical power 
of 6 X 10^^ erg s~^. Fig. [6] indicates that the jet forms a 
cylinder from the NE inner lobe to the northern middle lobe. 
If it is assumed that the power is uniformly distributed, then 
it is possible to estimate the amount of mechanical power 
deposited into a single grain and hence deduce the dust grain 
temperature using Eqn. @. 

We approximate the jet as a cylinder ^4' in diame- 
ter, which at the distance of Cen A is equivalent to ^4.2 
kpc. The cross-section of the cylinder at the distance of the 
northern dust cloud (^15 kpc) as seen from the nucleus of 
Cen A subtends an angle of ~16°. Hence the jet subtends a 
solid angle of ^0.005 sr at the distance of the northern dust 
cloud. This gives the jet a flux density of ^1.2 x lO^^erg s~^ 
sr~^. This equates to ^-^3.6 xlO~^'^W per grain. Balancing 
this power with Eqn (|4]) yields a dust temperature of ^2600 
K. If all the mechanical energy were dumped into the dust, 
the dust would not survive, so the transfer of jet mechanical 
energy would have to be incredibly inefficient. The southern 
cloud also poses a problem for this mechanism, since there 
is no evidence of a jet-cloud interaction in the region sur- 
rounding the southern cloud. Despite the striking alignment 
between the dust clouds and the jet, the jet appears to have 
little influence on the dust, suggesting that the alignment is 
purely a projection effect. 

We come to the same conclusion as Stickel et al. (2004). 
The simplest scenario that satisfactorily explains the origin 
of the dust is one in which the dust and gas in the clouds 
were once part of the late-type galaxy that merged with 
Cen A to produce its prominent dust lane. The only heating 
source that can account for the observed temperatures in 
both of the clouds is heating from the starlight of the evolved 
stellar population within Cen A. 



5 CONCLUSIONS 

The results of these HerschelSPlRE observations of Cen 
A have revealed two dust clouds external to Cen A and co- 
aligned with the axis of the AGN jet. The northern cloud lies 



© 2007 RAS, MNRAS 000, [T]-?? 



Herschel observations of Cen A dust clouds 15 



at a projected distance of ^15 kpc from the centre of Cen A 
and has been previously detected, while the southern cloud, 
at a distance of ^12 kpc, is a new detection in the sub-mm. 
We have used SED analysis and a single-temperature, mod- 
ified blackbody spectrum to fit the Herschel-SPIRE data 
and reprocessed MIPS data. The resulting fits show that the 
two clouds are thermal in origin and have roughly the same 
dust temperature and mass: Tnorth = 12.61^ 2 ^5 T^south = 
15.111:^ K; log(Mnorth/M0) = 5.8t°-j, log(Msouth/M0) = 
5.61q'2- The measured values for the northern cloud are 
consistent with previous measurements based on ISO data 
(Stickel et al, 2004). 

The clouds reside at the termini of the partial Hi ring 
that surrounds Cen A where the H i column density is high- 
est. The measured gas-dust ratios and Hi / H2 ratios are 
typical of a late-type gas-rich galaxy, which support the hy- 
pothesis that the dust and gas originated in the galaxy which 
has since merged with Cen A to form its dusty disk. Dust 
was not unambiguously detected in other parts of the Hi 
ring and we attribute this to the dust having a column den- 
sity below the SPIRE detection threshold rather than having 
a lower temperature or an absence of dust altogether. 

The dust clouds have also been detected in Hi, CO, 24 
\im and the southern cloud, at least, also exhibits 8 ^tm PAH 
emission. SFRs have been calculated based on UV-FIR trac- 
ers but the SFRs are too low to have produced the inferred 
mass of dust in a reasonable time frame. The low SFRs and 
the lack of UV emission in the southern cloud imply that the 
dust emission is unrelated to star formation. Dust heating 
by the evolved stellar population is the simplest mechanism 
capable of explaining the observed temperatures for both 
clouds. The jet therefore is unlikely to play a major role 
in the heating process and its alignment with the clouds is 
purely coincidental. 
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